Adiabatic quantum pumping and rectification effects in interacting quantum dots 
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We derive a formula describing the adiabatically pumped charge through an interacting quantum 
dot within the scattering matrix and Green's function approach. We show that when the tunneling 
rates between the leads and the dot are varied adiabatically in time, both in modulus and phase, the 
current induced in the dot consists of two terms, the pumping current and a rectification-like term. 
The last contribution arises from the time-derivative of the tunneling phase and can have even or 
odd parity with respect to the pumping phase cp. The rectification-like term is also discussed in 
relation to some recent experiments in quantum-dots. 



INTRODUCTION 



The idea of quantum pumping, i.e. of producing a dc 
current at zero bias voltage by time periodic modula- 
tion of two system parameters, dates back to the work 
of Thoulessi. If the parameters change slowly as com- 
pared to all internal time scales of the system, the pump- 
ing is adiabatic, and the average charge per period does 
not depend on the detailed time dependence of the pa- 
rameters. Using the concept of emissivity proposed by 
Biittiker et al4, Brouwer— related the charge pumped in 
a period to the derivatives of the instantaneous scatter- 
ing matrix of the conductor with respect to the time- 
varying parameters. Since then, a general framework to 
compute the pumped charge through a conductor has 
been developed for noninteracting electrons^. The inter- 
est in the pumping phenomenon has shifted then to the 
experimental^ investigations of confined nanostructures, 
as quantum dots, where the realization of the periodic 
time-dependent potential can be achieved by modulating 
gate voltages applied to the structure^. In case of inter- 
acting electrons the computation of the pumped charge 
becomes rather involved and few works have addressed 
this issue for different systems^ and in specific regimes. 
As for the case of interacting quantum dots, the pumped 
charge in a period was calculated by Aonc4 by exploit- 
ing the zero-temperature mapping of the Kondo prob- 
lem. A very general formalism was developed in Rcf. 9] 
where an adiabatic expansion of the self-energy based on 
the average-time approximation was used to calculate the 
dot Green's function while a linear response scheme was 
employed in Ref.[l4j]. More recently, another interesting 
study 15 was performed aiming at generalizing Brouwer's 
formula for interacting systems to include inelastic scat- 
tering events. 

In this work we present a general expression for the 
adiabatic pumping current in the interacting quantum 
dot in terms of instantaneous properties of the system 
at equilibrium, generalizing the scattering approach for 
noninteracting particles and discuss the limit of its valid- 
ity. To get a pumped current the two model parameters 
which are varied in time are the tunneling rates between 
the noninteracting leads and the quantum dot. In par- 



ticular, we let them vary both in modulus and phase 
through the adiabatic and periodic modulation of two ex- 
ternal parameters (e.g. gate voltages or magnetic fields) 
and show that a rectification-like term arises in the cur- 
rent due to the time-dependent tunneling phase. 

The plan of the paper is the following. In Sec[H]we in- 
troduce the model and relevant parameters. We develop 
the scattering matrix approach together with the Green's 
function formalism to derive the formula of the pumped 
current through an interacting multilevel quantum dot in 
the adiabatic regime at very low-temperatures. In Sec lIIII 
we specialize on a single-level quantum dot and give the 
explicit expression of the pumped current. Conclusions 
are given in Sec lIVI 



II. 



THE MODEL AND FORMALISM 



We consider a multi-level quantum dot (QD) coupled 
to two noninteracting leads, with the external leads being 
in thermal equilibrium. The Hamiltonian of the system 
is given by: 



H 



H. 



dot 



where H, 



J 



leads 
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with 



(cfc,<r,/3) the creation (annihilation) operator of 
an electron with spin a =],[ in the lead (3 — L,R 
and dispersion ep(k). The QD is described by the 



Hamiltonian H, 



dot 



Utia 



where 



djpdja with dj a (dj cr ) the creation (annihilation) 
operator of the electron with spin a and tj a the dot 
j-th energy level. The on-site energy U describes the 
Coulomb interaction. The tunneling Hamiltonian is 

given by H tun = Y,k,<r,i3j[ V K<rM( t ) c t,*,p d 3<> + H - c -h 
with time-dependent tunnel matrix elements Vk, a ,i3,j(t). 
For simplicity we assume that Vk,a,[3.j are spin indepen- 
dent, i.e. Vk,a,p,j = Vji,j and that both the modulus 
and the phase of Vpj(t) vary in time with frequency lu, 
i.e. Vpj(t) — \V/3j(t)\exp{i$p(t)}. Their explicit time 
dependence is determined by two external parameters 
(e.g. two gate voltages applied at the barriers of the dot 
or a gate voltage and a magnetic field) which are varied 
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adiabatically and periodically in time or by the presence 
of parasitic bias voltages. Two specific examples will be 
considered below. In particular, we will specialize on the 
case in which the tunneling phase can vary harmonically 
or linearly in time. The instantaneous strength of the 
coupling to the leads is instead characterized by the 
parameters T^ m (t,t) = 2irpVp, m (i)V^ n (t), where p is 
the density of states in the leads at the Fermi level. By 
varying in time VL,j(t) and Vnjfy) and keeping them out 
of phase, the charge Q pumped in a period T is related 
to the time dependent current I L (r) flowing through the 

left barrier, i.e. Q = J Q T (1tI l (t). 

While the exact formula for the current depends on 
time-dependent Green's function out of equilibrium, in 
the following we consider the adiabatic limit where the 
current depends only on the instantaneous equilibrium 
properties of the dot, i.e. on the retarded dot Green's 
function (GF). This situation is realized in the two fol- 
lowing cases. First, let us consider that only the mod- 
ulus of the tunneling matrix elements is varied in time, 
i.e. Vpj{t) = \Vpj(t)\exp{i&p}. Under the adiabatic 
condition, the tunneling rate varies slowly in time, and 
the quantum dot can be considered time by time in equi- 
librium with the external leads. The effect of quantum 
pumping is well described by an adiabatic expansion of 
the self-energy based on the average-time approximation 
as described in Ref . 0] and using the equilibrium relations 
to write the pumped current in terms of the retarded 
GF only. Let us consider now the situation in which the 
modulus of the tunneling terms is fixed while their phases 
are modulated in time, i.e. Vpj(t) — \Vpj\exp{i$p(t)}. 
This situation is equivalent (by a gauge transformation) 
to having a system biased with an ac external signal. 
In particular, the ac voltage applied to the leads is pro- 
portional to the time derivative of the tunneling phase 
dt$p(t) 19 . Since the tunneling terms are assumed to vary 
in time with frequency w, the ac signal forcing the quan- 
tum dot is proportional to the pumping frequency u> and 
thus can be considered as a small perturbation under the 
adiabatic condition. In particular, if one consider the 
case in which the tunneling phases vary linearly in time, 
*&/3(t) = iVt/2, this situation corresponds to an inter- 
acting quantum dot biased by a dc voltage V . Following 
the work by Meir and WingreeniS, the current / flow- 
ing through the interacting multilevel dot biased by a dc 
voltage V can be written as: 

I = (e/h) J de[f L (e) - f R {e)]Tr{G a T R G r T L K}, (2) 

where /l,r are the Fermi functions, Tl^r are the dot- 
leads coupling strengths and 7Z — Eq 1 S is the ratio 
between the fully interacting self-energy and the non- 
interacting one, responsible for the deviation from the 
Landauer-Biittiker formula (see Ref. [13], Eq.(lO)). In 
the zero-temperature limit and for a weak bias, E = So 
at the Fermi level and thus Eq. can written as 

I oc VTr{G a T R G r T L }, which corresponds to the usual 
linear response form, even though the Green's function 



are interacting ones. This argument is extensively dis- 
cussed in Refs.[l| (see Eq.s (37) and (38)) and [Icj]. 

Thus in general we expect that, when both the mod- 
ulus and the phase of the pumping parameters are var- 
ied in time, apart the usual dc pumping current a new 
term arises (that we call of rectification) which is pro- 
portional to the time derivative of the tunneling phase: 
I r = ^Tr{G a T R G r T L }d t ^ L - 

Since, as explained above, in adiabatic regime the cur- 
rent is determined by the instantaneous properties of the 
dot (retarded Green's function) and since in the zero tem- 
perature limit the usual linear response formula can be 
adopted for the calculation of the current, the pumping 
and rectification currents through the interacting quan- 
tum dot can be calculated by the scattering matrix ap- 
proach as well. In fact, as well known, the retarded 
Green's function is related to scattering matrix by the 
Fisher-Lee relation. We thus employ the scattering ma- 
trix formalism developed in Ref.[ll| where the charge 
current originated by an adiabatic pump is related to 
an expansion of the quantity {S(E,r)[f(E + id T /2) — 
f(E)]S'{E,T)} aa with respect to the time derivative op- 
erator id T /2 (here S(E,t) is the scattering matrix and 
f(E) is the Fermi function). The first order of this ex- 
pansion reproduces the famous Brouwer's formula^. Let 
us only stress that the scattering matrix formalism is well 
defined, not only in the noninteractin g ca se, but also for 
the interacting problem (e.g. see Ref.[l2|). The expres- 
sion of the pumped current (E, r) in terms of the time- 
dependent scattering matri x 4 ' 11 is: 

I P (E,t) = ^m{^S; a (E,r)9 T ^(£,r)}, 

a 

I p (r)= J ^ { -f)I^ E ,r), (3) 

where f(E) is the Fermi function and S(E, r) is the 
instantaneous S- matrix of the QD. It is given by the 
Wigner transform S(E, r) = J dte lEt S{r + t/2, r - t/2), 
where 

S a0 (t, = 5 aP 5(t - t') - Tr{M al3 (t, t')G r (t, t')}. (4) 

Here G r is the full retarded QD Green's func- 
tion, G£ >rn (t,f) = -iO(t - t')({d n (t),dl(t')}) and 
[M a/3 (t,t')] mn = 2mpV* trn (t)V(3, n (t'), In the limit of the 
pumping frequency w < T, i.e. under the adiabatic con- 
dition, the scattering matrix S(E, r) is expressed by the 
instantaneous Green's function of the dot as^: 

S a p(E,r) = 5 aP - 2nipJ2v:, m (T)G r nm (E,T)Vp, n (T). 

(5) 

When substituting ((SJ) into (J3J) to compute the current 
we need the time-derivative of the QD Green's function 
which satisfies the relation: 

d T G r {E, t) = G r (E, T)t r (E, r)G r (E, r), (6) 

where the dot symbol indicates a time-derivative and the 
matrix notation for the dot Green's function has been 
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used. The final expression obtained for the I t3 (E,r) for 
a multi-level quantum dot isi£: 

I^(E,t) = ^[Im{-iTr{t G r }}+Im{i($p-<f> ) x 
x Tr{F 3 G ( T' 5 G r }}+A /3 (£;,T)], (7) 

where 

A l3 (E, t) = 7m|-iTr{r /3 G a S r G r }} 

+ /m{rr{f /5 G r }Tr{r /3 G a }} 



x arMI^-ll^l)}- 



■ps 



(8) 



This expression has been obtained by considering explic- 
itly the time dependence of the modulus and phase of the 
tunnel matrix elements, and consequently of the leads- 
dot coupling function T(t). The symbol (3 stands for the 
L,R lead in correspondence of /3=R,L. The total dc cur- 
rent through the lead (3 is given by: 



2tt 



dE 



2-k/u 



drI^E,r)[-d B f(E)}. (9) 



The expression (|7J) represents our main result. It is valid 
for a multi-level QD and for any interaction strength in 
the zero temperature limit under the adiabatic condition. 
The first term in Eq.([7]) represents the pumping current, 
while the second one, proportional to the time-derivative 
of the tunneling rate phase, is the effective rectification 
term we have discussed above. It can also be written 
as Q(T)V e ff(r), where Q{t) cx Tr{T f) G a T^G r ] is the 
conductance of the structure, while V e ff(r) oc ($ ( g — $^). 
The last term in (J7J) contains information on the time 
derivative of the retarded self-energy and is zero for a 
single- level quantum dot within the wide band limit. 



III. TOTAL CURRENT FORMULA FOR A 
SINGLE LEVEL QD 

Up to now we have developed a theory of the pumped 
current valid in the case of a multi- level QD. We now 
specialize Eq. ([7]) to the case of a single level QD. Elim- 
inating the trace in ([7]) and considering the remaining 
quantities as c-numbers, the expression for the current 
simplifies to: 

I p (E,t) = -^-[i?e{f /3 G r } + r /3 |G r | 2 a T i?e{S r } + 

+(4 / 3-$ / 3)r /3 r^|G"'| 2 ]. (io) 

When the time-derivative of the tunneling phase is 
neglected the above formula is equivalent to the 
pumped current calculated by the self-energy adiabatic 
expansion 9 . 



In the following we consider the case of a single level 
QD both in the strongly interacting and non-interacting 
case and describe the behavior of the charge Q (in unit 
of the electron charge q) pumped per cycle in the zero- 
temperature limit. The Fermi energy /i is set to zero as 
reference energy level, while the static linewidth T ' is 
assumed as energy unit (typical value for Tq is 10 fj,eV). 
In the noninteracting case, i.e. when the QD Green's 
function becomes a scalar, the expression for the instan- 
taneous pumping current is explicitly given by: 



2tt 



(E - e )r? 



{E - e ) 2 + (T/2) 2 



(E-e r + (T/2y 



(11) 



When <E>^ — $^ = 0, the charge pumped is zero when the 
level is resonant (e = 0). 

In the case of a strongly interacting quantum dot, i.e. 
in the infinite-?/ limit, we take the expression of the QD 
Green's function as in Ref.|13l. The current is: 



Jl_ 

2tt 



(E - Eo )(l - n 9 )tP 



(£-e ) 2 + (l72)2(l-n ff )2 
r^($ / 3-4> /3 )(l-n ff ) 2 



(£-£ ) 2 + (r/2)2(l-n ff )2 



(12) 



where the occupation number n a on the dot has to 
be determined self-consistently by the relation n a — 
(2m)- 1 J dEG< a {E), where G< a is the QD lesser Green's 
function. 

In order to show the effects of the time-dependent tun- 
neling phase we report below the numerical results of the 
charge Q pumped per cycle. The pumping cycle is de- 
termined by the periodic time variation of the leads-dot 
coupling strength, where T a (t) = Tg + V" s'm(uit + (p a ), 
while for the tunneling phase two cases can be considered. 
Either it varies harmonically $ Q (t) = $^ sm(ojt + ip a ) 
with the same frequency of the two external gate volt- 
ages (this case is shown in Fig[TJ) or it varies linearly in 
time <&Q.(t) = 3>£t, e.g. when a parasitic gate voltage is 
present, (this case is shown in FigEJ) . The quantity ip a 
is the pumping phase that we take different from zero 
between L and R lead. 

In FigQ] we plot Q as a function of the energy level Eq 
by fixing the other parameters as: ji = 0, <p = tt/2, 



t* = 0.2, r£ = o.4, r£ = r# = i,n = 0.3, = 0.2. 



In particular, in the upper panel, the charge induced by 
the pumping (triangle) and the charge due to the rec- 
tification term (box) is shown for a non-interacting dot 
(U = 0). The total charge (empty circle) is significantly 
modified by the presence of the rectification term which 
is a non- vanishing quantity at the Fermi energy (eo = 0). 
The lower panel in FigQ] shows the behavior of the rec- 
tification current as a function of Sq in the strongly in- 
teracting limit (U — > oo) and by choosing the remaining 
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parameters as in the upper panel. While the general as- 
pect of the total charge (empty circle) is only marginally 
modified by the strong correlations in the specified region 
of parameters, the rectified charge (box) shows a pro- 
nounced asymmetric behavior with respect to the level 
of the dot. In Fig|2] we focus on the case of time- linear 
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FIG. 1: Charge Q pumped per pumping cycle as a function 
of the dot level eo. The pumping (triangle) and the rectifica- 
tion (box) contribution to the total charge (empty circle) are 
shown for the non interacting case (17 = 0) in the upper panel 
and for strongly interacting dot (17 — > oo) in the lower panel. 
Both figures are computed for the following choice of param- 
eters: n = o, ip = tt/2, it} = 0.2, r£ = 0.4, r£ = r£ = i, 

= 0.3, = 0.2. The pumping cycle is determined by: 
F a (t) = rg 1 + rZ sin(oj/j + ip a ),$ a (t) = sin(wi + (p a ), with 
a — L,R and <pl = 0, tpa = <p 

variation of the phase $ Q (i) = and take the param- 
eters as in Fig[T] In the upper panel, the pumped charge 
Q is plotted as a function of the dot level eo m the non- 
interacting case (U = 0). Contrary to the previous case, 
the rectification contribution (box) is dominant over the 
one induced by the pumping mechanism (triangle) and 
thus the total charge (empty circle) is mainly affected by 
a resonant-like behavior. In the lower panel, the results 
for U — ► oo are shown. Apart from a renormalization 
of the linewidth of the resonance induced by the factor 
(1 — riff) in the numerator of Eq. (|12[) . a behavior similar 
to the one of the non- interacting system is found. A slave 
boson treatment with the inclusion of a renormalization 
of the dot energy level, could in principle modify this pic- 
ture. Let us note that when the tunneling phase varies 



harmonically both the pumping current and the recti- 
fication current follow the same sin(yj) behavior w.r.t. 
the pumping phase (p. In Fig[3]we show the behavior of 
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FIG. 2: Charge Q pumped per pumping cycle as a function 
of the dot level eo in the non interacting case (upper panel) 
and strongly interacting case (lower panel). The pumping 
(triangle) and the rectification (box) contribution to the total 
charge (empty circle) are shown in the zero temperature-limit 
by setting the remaining parameters as follows: fi = 0, <p = 

tt/2, rj = 0.2, rt = 0.4, r£ = r? = 1, = 0.3, = 0.2. 

Differently from FigfTJ the pumping cycle is determined by: 
r a (t) = I^ +r£sin(^ + tp a ), $ Q (i) = $£t, with a — L,R 
and ipL = 0, (fR — tp. 

the charge induced by the pumping term (upper panel) 
and by the rectification (lower panel) with respect to the 
pumping phase ip and by fixing the dot level to eo = 0.3 
and the remaining parameters as done in Figj2] The full 
line in both panels represents the result for the U = 
case, while the full circles (•) represent the curves com- 
puted for the infinite-?/ case. Let us note that while the 
pumping term follows the conventional sin(</?)-behavior 
as in Brouwer theory, the rectification contribution takes 
the form V e f f[A+ B cos(ip)], where the coefficients A and 
B for the U = case are explicitly given by: 



A 



-nL-nR. 
1 1 



B = - 
with £ = (/i — Co 



C + (Fo/2) 2 

rMrgxr^ - r«) 2 + 2t^{a^ - vl) 



(13) 



32[£+(r /2) 2 ]3 
2 . The different symmetry of the pump- 



5 



ing and rectification current has already been reported 
in experimental works in quantum dots2i and some theo- 
retical explanations have been proposed^ 2 -. Furthermore, 
the analysis of the coefficient B shows that the charge 
transferred by the rectification effect can be significantly 
increased by coupling the dot region to the leads in an 
asymmetric way (r^ ^ Tq), e.g. by using tunnel barriers 
with very different transparencies. The above results re- 
main almost unchanged in the strongly interacting case 

(tr->oo). 

In the limiting case in which the phase difference be- 




1 2 3 4 5 6 



-0.175 
-0.180 




1 2 3 4 5 6 



FIG. 3: The behavior of the pumping current (upper panel) 
and of the rectification current (lower panel) with respect to ip 
is reported for so = 0.3 and taking the remaining parameters 
as in Fig[2] Notice that the pumping term is proportional to 
sin(y>), while the rectification one is proportional to cos(<p). 
Both in the upper and lower panel, the full line represents the 
non interacting result, while the interacting case is indicated 
by full circles (•). 



tween the tunneling barriers is kept zero, i.e. <pi, = (pji — 
0, the pumping term is exactly zero and the rectification 
contribution acts as a quantum ratchet 26 . 



IV. CONCLUSIONS 

Within the Green's function and scattering matrix ap- 
proach we have analyzed the quantum pumping current 
through an interacting quantum dot when both the mod- 
ulus and the phase of the model time-dependent parame- 
ters, in our case the leads-dot tunneling rate, is adiabati- 
cally varied. In this way it has been possible to derive an 
expression for the pumped current containing an effective 
rectification term due to the time-dependent phase. Such 
contribution can be written in a Landauer-Buttiker-like 
form, even though for the interacting system, when the 
zero temperature limit and the adiabatic conditions are 
met. The numerical analysis also show that when the 
tunneling phase varies linearly in time the rectification 
term is even with respect to the pumping phase tp, i.e. of 
the form a + b cos(tp) , in contrast to the usual pumping 
term which is odd. The mentioned contribution could be 
related to the experimentally observed rectification ef- 
fects in quantum dots 2 ^. 

In particular, we have been considering an open system, 
but concerning closed systems (e.g. annular devices with 
a quantum dot), the tunneling phase contribution to the 
pumped charge could find a natural interpretation in a 
complex phase of geometric nature^ 5 . This phase would 
be a Berry phased. Thus the detection of a rectification 
current in addition to the pumping one could be an in- 
direct probe of a Berry phase. 

The proposed analysis could be easily generalized up to 
the second order in the pumping frequency uj allowing to 
describe features involved in the moderate non-adiabatic 
limit. 
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